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ABSTRACT 

We study global non-axisymmetric stationary perturbations of aligned and unaligned loga- 
rithmic spiral configurations in an axisymmetric composite differentially rotating disc system 
of scale-free stellar and isopedically magnetized gas discs coupled by gravity. The infinitely 
thin gas disc is threaded across by a vertical magnetic field with a constant dimension- 
less isopedic ratio A = 2'K\fGTS^'^ j of surface gas mass density S*^^) to B^ with G being 
the gravitational constant. Our exploration focuses on the relation between the isopedic ra- 
tio A and the dark matter amount represented by a gravitational potential ratio / = $/$ in 
order to sustain stationary perturbation configurations, where $ is the gravitational potential 
of a presumed axisymmetric halo of dark matter and $ is the gravitational potential of the 
composite disc matter For typical disc galaxies, we explore relevant parameter ranges nu- 
merically. High and low A values correspond to relatively weak and strong magnetic fields 
given the same gas surface mass density, respectively. The main goal of our model analysis is 
to reveal the relation between isopedic magnetic fields and dark matter halo in spiral galax- 
ies with globally stationary perturbation configurations. Our results show that for stationary 
perturbation configurations, fairly strong yet realistic magnetic fields require a considerably 
larger amount of dark matter in aligned and unaligned cases than weak or moderate magnetic 
field strengths. We discuss astrophysical and cosmological implications of our findings. For 
examples, patterns and pattern speeds of galaxies may change during the course of galactic 
evolution. Multiple-armed galaxies may be more numerous in the early Universe. Flocculent 
galaxies may represent the transitional phase of pattern variations in galaxies. 

Key words: galaxies: haloes — galaxies: ISM — galaxies: kinematics and dynamics — 
magnetic fields — magnetohydrodynamics (MHD) — waves 



1 INTRODUCTION 

Since the seminal studies of the galactic density wave theory in the 1960's pioneered by Lin & Shu (1964, 1966), many advanced works 
have been done to describe the large-scale dynamics of 'grand-design' spiral structures in disc galaxies. In terms of theoretical development, 
it is natural to first concentrate on linear perturbation theories. As synchrotron radio observations have revealed since the early 1970's that 
magnetic fields are ubiquitous in disc galaxies (e.g. Sofue, Fujimoto & Wielebinski 1986; Beck et al. 1996; Zweibel & Heiles 1997; Brown, 
Taylor, Wielebinski & Mueller 2003), magnetic fields have been recognized as an important aspect in the model treatment of spiral galaxies 
and thus the magnetohydrodynamic (MHD) density wave theory was developed (Fan & Lou 1996; Lou & Fan 1998). Over the past decade, 
a research area has been developed on the so-called scale-free discs and perturbation structures therein (e.g. Syer &l Tremaine 1996; Shu et 
al. 2000; Shen & Lou 2004; Lou & Zou 2004; Shen, Liu & Lou 2005; Lou & Zou 2006; Wu & Lou 2006; Lou & Bai 2006). The scale-free 
condition corresponds to a power-law dependence of all physical quantities on the cylindrical radius r, for example, the background surface 
mass density is Eo oc 7-~2/3-i p being a constant scaling index. Such analysis was first initiated for relatively simple models of only one 
gas disc. We are now able to explore a composite system of a stellar disc and an isopedically magnetized gas disc embedded in a dark matter 
halo. For this overall configuration, we study possible large-scale perturbation structures and stationary MHD density waves. To simplify the 
mathematical treatment, the stellar disc is approximated as a fluid and the magnetized gaseous disc as a magnetofluid and both are assumed to 
be geometrically razor-thin and scale-free discs. In general situations, the scaling exponent (3 may vary radially in a certain range for different 
disc galaxies. For more or less flat rotation curves which can be grossly determined in many spiral galaxies from large radii to the centre (e.g. 
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Rubin 1965; Roberts & Rots 1973; Krumm & Salpeter 1976, 1977; Bosma 1981; Kent 1986, 1987, 1988), index f3 may be regarded as almost 
zero. According to the classical Newtonian gravity theory, a significantly larger amount of mass than the visible mass of the disc is needed 
in order to sustain observed grossly flat rotation curves. As a consequence of this discrepancy between the visible and the required mass as 
well as the disc stability property, the presence of a massive dark matter halo has been proposed (e.g. Ostriker & Peebles 1973; Ostriker et 
al. 1974; Binney & Tremaine 1987). Nowadays, the inclusion of a massive axisymmetric dark matter halo is a common practice for studying 
large-scale galactic dynamics. 

Meanwhile, observations of disc galaxies nowadays can be used to estimate open magnetic field strengths and configurations on both 
sides of a galactic disc plane. While these open magnetic fields are most certainly interlaced with closed coplanar magnetic fields, at this 
stage of our theoretical model development, we consider only open magnetic field crossing the disc plane almost ver ticallyQ Shu & Li 
(1997) introduced the so-called isopedic condition of Bz{r, 9, t) oc E'®' where is the z— component of the magnetic field and E'®' is 
the gas disc surface mass density. The ratio E'*'/_Bz thus remains constant in time and space. This situation is referred to as the isopedic 
magnetic field configuration. In fact, this constant ratio can be actually inferred by the magnetic field-freezing effect according to the ideal 
magnetohydrodynamic (MHD) equations (Lou & Wu 2005). Effective modifications to the gravitational potential and the gas pressure are 
derived in Shu & Li (1997) for a single isopedically magnetized gas disc. This isopedic condition was first proposed as an ansatz based 
on numerical simulations of magnetized cloud core formation through ambipolar diffusions (e.g. Nakano 1979; Lizano & Shu 1989; and 
references cited in Shu & Li 1997). On the basis of earlier work of Shu & Li (1997), Lou & Wu (2005) generalize their results into theorems 
for two gravitationally coupled discs where only the gas disc is isopedically magnetized. Moreover, they actually demonstrate that a constant 
ratio A = 27r\/GE'^'/_B2, with G being the gravitational constant, can be proved as a consequence of frozen-in condition by the standard 
ideal nonlinear MHD equations. In this model analysis, we shall adopt the same expressions of Lou & Wu (2005) for the effect of an isopedic 
magnetic field anchored in the gas disc component. 

In order to catch the essence of MHD density wave equations from analytical and numerical explorations, a few assumptions and 
simplifications are invoked such as the stationary condition in our frame of reference. As first advanced and hypothesized by Lin & Shu 
(1966), the quasi-stationary spiral structure (QSSS) has been studied since the early 1960's (e.g. Bertin & Lin 1996 and extensive references 
therein). There the pattern speed of density wave travels very slowly as compared to the differential rotation speed ve of the disc matter. Our 
stationary requirement can be simply regarded as a limiting case of the QSSS approximation. Based on observations of the nearby spiral 
galaxy M81, Visser (1980a, b) has shown that the density wave theory does indeed account for gross aspects of the large-scale dynamics of 
this spiral galaxy. 

Kendall et al. (2008) presented further observational data analysis of the nearby galaxy M81, where the large-scale spatial offset between 
the (old) stellar spiral wave and the gas shock front was examined. Theoretical considerations (e.g. Gittins & Clarke 2004; Chakrabarti 2008) 
have shown that the offset between the gas shock front (i.e. the maximum gas density) and the spiral potential minimum (determined by the 
spiral pattern in the old stellar disc component) may yield clues about the lifetime of a spiral pattern. By comparing observational results 
with the predictions of Gittins & Clarke (2004), Kendall et al. (2008) found that the galaxy M81 probably possesses a long-lived spiral 
pattern. Very recently, Chakrabarti (2008) performed numerical simulations to demonstrate that an angular momentum transport takes place 
from gas to stars in a composite system, leading to an approximately time-steady spiral structure in the stellar disc component (treated as 
coUisionless N-body particles). A major part that plays an important role is the relatively cold gas disc particles modelled by smoothed 
particle hydrodynamic (SPH) and is dissipative due to an artificial viscosity. Regarding the spiral structure in the gaseous disc component, 
Chakrabarti et al. (2003) have shown, that the so-called ultraharmonic resonances that the gas experiences in a fixed stellar spiral potential 
can lead to highly time-dependent gas responses, not allowing long-lived spiral patterns in the gas disc. In our model of a composite disc 
system, there is no fixed stellar potential but rather stellar and gaseous discs are gravitationally coupled to each other in terms of dynamical 
adjustment. While the simulation results reveal only a long-lived stellar spiral pattern, for our semi-analytical calculations, a stationary 
requirement for both spiral patterns in stellar and gaseous discs may be imposed without inconsistency. This offers a useful construct for 
modelling large-scale patterns in isopedically magnetized spiral galaxies. 

Theoretically, two classes of stationary MHD density wave solutions with in-phase and out-of-phase density perturbations are expected 
to exist in a composite disc system according to Lou Sl Shen (2003). The physical reason for these two different couplings between the stellar 
and gaseous discs is a dynamic interplay between star formation and gas clumping processes (Shen & Lou 2004). When the surface mass 
densities are out of phase, the self-gravity of the composite disc system is reduced. The azimuthal propagation of density wave is therefore 
faster and if a stationary configuration is desired, then the disc should rotate faster. When the surface mass densities are in phase, the self- 
gravity is enhanced. The azimuthal propagation of the density wave is slower and the disc should then rotate slower in order to maintain a 
stationary configuration (Lou & Shen 2003). 

Our present investigation mainly focuses on the variation of the amount of dark matter in the halo needed to sustain stationary perturba- 
tion configurations as the isopedic magnetic field strength varies. Therefore, we use typical values of disc galaxies for different parameters 
such as the 'sound speed' or the surface mass density. By invoking several physical requirements for the dark matter halo such as a non- 
negative potential ratio and so forth, the physical parameter range of stationary perturbation solutions can be identified. 

In Section |2l the problem formulation, the basic formulae and the coupled fluid and magnetofluid equations for the two discs are 
presented. In Section [3] we discuss effective modifications to the gravitational potential and the enthalpy in the gas disc when an isopedic 
magnetic field is introduced. In Section |4l we derive several relevant relations for the background rotational equilibrium configuration. In 
Section[5] MHD perturbations are introduced and the relevant global stationary dispersion relation is derived. For globally stationary MHD 
perturbation configurations, we calculate the dark matter amount and discuss which parameters affect this amount. Section|6]describes briefly 
the values for the parameters adopted for this model analysis. In Section|7l we study aligned cases where perturbations of different radii can be 



The studies of coplanar magnetic fields in a rotating thin disc system can be found in Lou & Zou (2004, 2006) and Lou & Bai (2006). 
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aligned radially. In Section[8l we study unaligned logarithmic spiral cases with a constant radial flux of angular momentum. Such unaligned 
cases have MHD perturbations which are systematically phase shifted azimuthally as a function of radius r. In unaligned cases, MHD density 
waves propagate in both azimuthal and radial directions. The constant radial flux of angular momentum is realized by setting the perturbation 
scale-free index /3i to 1/4 (e.g. Goldreich & Tremaine 1979; Shen et al. 2005). In Section|9] we discuss two limiting cases of no magnetic 
field and only a single isopedically magnetized gaseous disc. In the last Section 10, we summarize our work and draw conclusions. 



2 NOMENCLATURE AND BASIC FLUID AND MAGNETOHYDRODYNAMIC (MHD) EQUATIONS 

We first introduce notations for the relevant physical variables used in our model formulation and analysis. We largely adopt the same 
nomenclature of Wu & Lou (2005) for a model of a single scale-free disc with an isopedic magnetic field. Whenever the superscript (i) is 
attached to a physical variable, the relevant equation is valid for both stellar (s) and gaseous (g) discs in parallel. Cylindrical coordinates 
(r, 8, z) are adopted for our model formulation and we confine our analysis within the z — plane for an infinitely thin composite 
disc system in differential rotation. For example, variable E^'''(r, 6, t) is the two-dimensional surface mass density of the stellar disc and 
S'®' (r, 9, t) is that of the isopedically magnetized gas disc. The total surface mass density E(r-, 9, t) = E^"^ + E'"' is the sum of these two 
disc surface mass densities. Two different coplanar bulk flow velocities are possible in the two discs: the radial bulk flow velocity Ur'' (r, 9, t) 
and the azimuthal bulk flow velocity Vg^^ (r, 6, t). For the background rotational equilibrium disc configuration of axisymmetry, the radial 
flow velocities should vanish. The two-dimensional pressure (i.e. vertically integrated pressure) in each disc is described by n''-'(r, 9, t). 
Between the two-dimensional variables and their three-dimensional counterparts, the following integral relations exist, 

n^'\r, 9, t) = J p^'\r, 9, z, t)dz and E<''(r, 9, t) = j p^'\r, 9, z, t)dz , 

Az Az 

where p*^*' is the three-dimensional pressure and p^-^^ is the three-dimensional mass density, and Az is the disc thickness presumed to be 
very small compared to radius r. Thus, II'*' is the vertical integral of pressure p'*'. The three-dimensional polytropic equation of state is 
p*-'' — k^'' [p*^*']^ with 7 being the polytropic index and k^^ being related to the disc entropy. The vertically integrated barotropic equation 
of state can then be expressed as 



(1) 



(2) 



where index n > is the barotropic index, fc'*' > for a warm disc and fc*^'' = for a cold disc. The barotropic sound speed a(') in each 
disc is then defined by 

p*'] ^ = dn^') /dE^'' = nfc''^ [e^''] = nn^*' /E^'' . (3) 
The 'sound speed' a^"' in the stellar disc is effectively related to the velocity dispersion of stars. The enthalpy Z/*^*' in each disc is given by 

J E(>) (n - 1) (n - 1) ^ ' 

and is proportional to the square of 'sound speed' in each disc. The gravitational coupling between the stellar and magnetized gaseous discs 
is described by the well-known Poisson integral 

where G = 6.67 x 10^* g^^ cm'' s^^ is the gravitational constant. The gravitational potentials of the two discs can be separately written as 

with $ = + In a disc galaxy, an axisymmetric massive dark matter halo is represented by a gravitational potential $. A recent 
numerical exploration (e.g. Diemand et al. 2008) has shown that the dark matter halo is not smooth and uniform, but rather consists of many 
smaller clumpy subhalos. For simplicity, we assume the dark matter halo to be grossly axisymmetric on large scales and unperturbed in our 
model analysis at this stage. The composite system of two coupled discs with an axisymmetric dark matter halo is described by coupled 
nonlinear fluid equations Q — (|9| below. 

9E(') ^ 1 a r ^r,. r,.l . 1 d 



dt r dr 



dt ^ dr r 89 r E(*) dr dr dr ' 

dv^ ^dv^ ve^_ 1 + $ + ¥] 

dt dr ^ r d9 ~^ r ^ E(') rd9 rd9 ^ rd9 

(e.g. Shen & Lou 2004). Equation Q is the mass conservation, and equations {Sj and ^ are the radial and azimuthal momentum equations. 
These fluid equations for both stellar and gaseous discs are confined to the two-dimensional disc plane at z = 0, while Poisson integral 
equation ([5} is three-dimensional. These equations do not involve an isopedic magnetic field at this stage. 



4 M. Xiang-Gruess, Y.-Q. Lou, W. J. Duschl 



2.1 Scale-Free Discs and Relevant Parameter Ranges 

Scale-free discs are represented by power-law forms with a few disc index parameters. We write a self-consistent form of disc solution as 



vi"'^ ^ e^^\ip)r~^ , 4"' = e'^'Mr"" , (10) 

vl"^ ^b^^\i^)r-'^ , u^"' = &<"^(v9)r-^ , (11) 

E(«' = S'-'\^)r'^^-^ , = S^-'\^)r-^^-^ , (12) 

n{9) , n'^' =fc(^)[E(«)]", (13) 

H^'^ = Q'^'\^)r-^^ , H^'^ = Q^'\ip)r-^'' (14) 

<!> = -P((^)r"^'' , ¥ = -Pr-^'^ , (15) 



where e'^', fe'"', S'^'\ Q<^', and P are functions of only ip which is an argument abbreviation ip = 6 + /ilnr with /i being a parameter, and 
P for the dark matter halo is a constant coefficient independent of ip. The detailed mathematical procedure of constructing such self-similar 
scale-free solutions can be found in Lynden-Bell & Lemos (1999). By comparing equations Q and il4\ . we immediately obtain a relation 
between the two power-law indices n and /3, namely 

(2/3 + l)(n- 1) =2/3, ^ n = (l+4/3)/(l + 2/3). (16) 

Based on expressions ^ and l|3j, we find that for warm discs with fe'*' > 0, the barotropic index n must also be positive in order to ensure a 
positive right-hand side (RHS) of equation l[3j (i.e. a real sound speed). This physical requirement for warm discs leads to either (3 > —1/4 
or /3 < — 1/2. For cold discs with fc*-'^ ^ 0, there is no constraint on the barotropic index n. Another empirical constraint which arises from 
observational results is that E*^*' should decrease with increasing r which requires a (3 larger than —1/2. With the assumption of decreasing 
E'^', a singularity of mass would arise in the central region as r — > 0^. This problem can be remedied by requiring a finite integral 

^1-2/3 



and this then leads to /3 < 1/2. Therefore for warm discs with fe''' > 0, the scaling index range is /3 G (—1/4, 1/2) and for cold discs with 
A:''' = 0, we simply require /3 € (-1/2, 1/2). 



f-r PZ7T fZTT 

lim / / E'*'(r,e»)rdrd6l = / S''*>(6i)de» lim 
>— »ojo Jo Jo *" 



3 AN ISOPEDIC MAGNETIC FIELD ACROSS A SCALE-FREE THIN GASEOUS DISC 

Effective modifications of the gravitational potential and the pressure in a composite system of two discs due to the very presence of an 
isopedic magnetic field were derived in reference to the study of singular isothermal discs (SIDs) by Lou & Wu (2005). SIDs represent only 
a special class (i.e. /9 = 0) of scale-free discs in our more general formalism here. We now briefly summarize their results and introduce the 
relevant variables in the current model context. Following Li & Shu (1996), we define a constant dimensionless ratio A as 

A = 2ttG^^'^E'-^'> /B^ = = constant > . (18) 

Parameters A and A indicate how strong the magnetic field is as compared to the surface mass density of the gas disc; they are proportional 
to B^^ such that for weak or strong magnetic fields they become very large or small, respectively. Parameter 77 is a dimensionless ratio of 
the total horizontal gravitational acceleration |/||| = + /y^'l (continuous across the thin disc along the vertical direction) to the total 
vertical gravitational acceleration just above the two discs with a total surface mass density E = E<»' + E<^\ namely 

= \hL ^ (19) 

' ~ 27rGE 27rGE(9) ' 

where the gaseous disc contribution to /|| is f:^^ — — V||<E>'^'. The two-dimensional gradient operation Vy within the disc plane coincident 
with 2 = and in terms of cylindrical coordinates (r, z) is simply 

^'^-^-^a^<^ ^ (20) 

or r oO 

where eV and eo are unit vectors along the radial and azimuthal directions, respectively. The sum of the magnetic tension force and the 
horizontal gravity force acting on the isopedically magnetized gaseous disc is 

= fin + .ft ^ = ^./]^' , (21) 

with as the magnetic tension force acting in the magnetized gaseous disc and e = 1 — 1/A^. In other words, the magnetic tension force 
fin and the horizontal gravity of gaseous disc /^^^ act always in opposite directions. Effectively, equation M\\ leads to the modification of 
the gravitational potential for the gaseous disc as 

${9) ^ , (22) 
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where e is referred to as the reduction factor because it is always less than unity and may become negative. The situation of e < happens 
when the magnetic tension force overwhelms the gas disc horizontal gravity. Due to the additional magnetic pressure, the total pressure in 
the magnetized gaseous disc is effectively enhanced by 

1 -I- ^ 

n(9)_^en(9\ where 9 = 1+ '' >0, A=(1 + 5"')A, 5 = e'^'/S^'^ (23) 

Here, Q is referred to as the enhancement factor because we always have 6 > 1. For the effective modification of disc enthalpy H^-^\ we 
assume for simplicity that O is independent of the cylindrical coordinates r and 9 and then get iif'®' — > OH^'^K 

In summary, the gravitational potential cl?'^', the pressure n*^®' and the enthalpy //'^^ in an isopedically magnetized gaseous disc are 
effectively modified according to — > e<l?'®\ n'*' BIl'®' and H^^^ — » while the counterparts of these variables remain 

unchanged in the stellar disc. The magnetosonic speed squared in the isopedically magnetized gaseous disc is then 



= „efc(s) ^e'^'J = e p»>J . (24) 

In reference to equations (|7]) — ([9]), an isopedic magnetic field in the gaseous disc leads to the following modified set of coupled equations. 



9S« ^ 1 a_ r ^ti^ ml . 1 d 
dt r dr 



H"V'']+;^[S«.«]=0, (25) 



dt ^ dr r 89 r Ef'') dr dr dr ' 

dt dr r d9 r E(s) dr dr dr ' 

"^T"^^" ~&r^~ r ^~Ew7~ae 7d9 ^ 7d9 ' ^^^^ 

dt ~^ dr ~^ ~ d9~ r ^ ~S(s)r d9 Vd9 ^ Vd9 ' ^^^^ 

Equations (|25j — (1291) bear very similar form of hydrodynamic equations Q — © with the effect of an isopedic magnetic field being 
subsumed into two dimensionless parameters B and e related to the magnetized gaseous disc component (Shu & Li 1997; Wu & Lou 2006). 



4 EQUILIBRIUM CONFIGURATION OF A COMPOSITE ROTATING DISC SYSTEM 

In a stationary equilibrium, the stellar and gas discs in a composite disc system rotate with different angular speeds in general while satisfying 
the basic nonlinear fluid-magnetofluid equations. For both discs in coupled rotational equilibrium of axisymmetry, the gravitational accel- 
eration caused by the dark matter halo and the two discs together is the same in the two radial force balance conditions [i.e. eqns ( I30t and 
l |31t below]. Meanwhile, the gas pressure and magnetic Lorentz forces together in the gaseous disc are different from the effective pressure 
force produced by stellar velocity dispersion in the stellar disc in general. These naturally lead to two different angular speeds of the two 
discs. Equilibrium variables of a composite rotating disc configuration are denoted by a subscript 0. For a rotating disc configuration in a 
stationary axisymmetric equilibrium, radial velocities vanish with v^.^ = v^'q = and Vgg^ = fig^V, Vg^g = Og^'V with fig*' = Q^^r) 
being the differential angular rotation speed of each disc. The mass conservation equations are consistently satisfied by our prescription. 
Radial momentum equations l |26t and i27\ then lead to two balance conditions 

_a[i?r+*^^^+-f^l+*ol, (30) 



r dr 



r dr 
In general, the radial derivative of disc enthalpy is related to the 'sound speed' as follows 



(31) 



dr dr 



(n-l)L-o 



(2/3 + 1) . (32) 



For the total surface mass density in a rotational equilibrium configuration, we simply take Eo = Eg^'^+Eg"-' = 5g^V~'^''~^+5'o''V~^''~^ — 
S'gr"^''"^ where Sg"' and S'g®' are two positive constant coefficients (see also equation [T2]l and So = Sg"^ + Sg"' is the sum of these two 
coefficients Sg^' and Slf \ Prescribed as such, one can make use of formulae in Qian (1992) for the total gravitational potential 

•l-o = -GrEoFo(/3) = ~GYo{(3)Sor-^'' , (33) 
where the coefficient factor Yo{P) is related to the standard F functions by 

y (n. _ nr{l/2-mP) .... 
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For later global non-axisymmetric perturbation analysis, we also introduce below a generalization of yo(/3) as Ym{P) defined by 

, . ^ ^r(m/2-/?+l/2)r(m/2+/3) 

- r(m/2 - /3 + l)r(m/2 + (3+ 1/2) ' ^ ' 

where m is an integer in the complex phase factor exp(— im&) for characterizing non-axisymmetric coplanar perturbations. As shown by 
Qian (1992), the valid range for expression l l35t is — m/2 < l3 < (m + l)/2. For the smallest m = 1, the valid range of /3 is (—1/2, 1) 
which is wider than the derived (3 range (—1/4, 1 /2) for warm background discs; and therefore we do not worry about the parameter regime 
of [3 in our subsequent global perturbation analysis. This leads to the following relation for the horizontal-to-perpendicular gravity ratio rj as 
defined by equation l |19t , viz. 

^ l-Vll-fol ^ - 2GpYo{l3)r-'^-'Soer\ ^ PMP) = , .... 
^ 27rGEo 27rGSor-2/3-i tt ~ ' ^ ' 

For each equilibrium configuration in the composite disc system, the gravitational potential and its first radial derivative are given by 

$« = -GrE«Fo(/?) = -GFo(/?)r-'''5'« and Mi = 2/30^0 (/?) S« . (37) 

To measure the effect of a dark matter halo, we introduce a gravitational potential ratio parameter / as 

f^ , ^° (38) 

for the ratio between the dark matter halo potential to the background potential of the two coupled discs together. The radial derivative of the 
axisymmetric dark matter halo potential can now be simply expressed as 

^ = — = .f2f3GYo{l3)Eo . (39) 

or or 

With the expressions for enthalpies and gravitational potentials derived above, radial force balances OOb and J31b become 

Ho]' + [«*^'] '(2/3 + 1) = 2l3rGYo{P) [E*^' + e'^'] (1 + /) , (40) 
[vil^] ' + e '(2/3 + 1) = 2PrGYo{f5) { [E^ + E'''^] (1 + /) - (1 - e)E('" } . (41) 

The physical properties of the two discs are related by equilibrium conditions ( I40t and ( I4U . The two polytropic sound speeds a'^' and a'^^ 
in the axisymmetric background stellar and gaseous discs are respectively 

_ „i.(s) l^(s)yi-l _ _^j^(s) /^(s) _ {I + 4/3) 2/3/(1+2/3) 

(1 + 2/3) 

[a(«'] ' = nfc(^' [e(«'] = nUlf^ /E^^' = ^i^fc'^' [si'^] ' . (43) 

(1 + zpj 

Combining equations < 14Qb and i HTt . one can readily deduce the equilibrium surface mass densities of two discs as 

^ [«^oT+[a'->r(2/3 + l) 

" 2/3rG11,(/3)(l + /)(l + 5o) ' ^ ^ 

^ [.<gT + e["'^'r(2/^ + i) 

° 2/3rGyo(/3)(/ + 1/So + //5o + > 



(45) 



with So = Eq^^/Eq"-* being the ratio of gaseous to stellar disc surface mass densities. As the surface mass densities are positive, we obtain a 
necessary inequality between the potential ratio / and the ratio A by using the fact that the quantities So, Yo{f3) and r are all positive 

/ + l/5o + f/So + e > for /3 > . (46) 

Note that e = 1 — 1/A^ may be negative for a strong isopedic magnetic field. This is the first requirement on potential ratio / parameter 
which will be used later. Since (3 > —1/4 for all discs, the numerators of equations J44t and ( 145 1 are always positive which leads to the 
condition that /? in the denominator must be positive since all other factors are positive. This condition restricts our exploration to the range 
of/3 = (0, 1/2). 

5 TWO-DIMENSIONAL COPLANAR MHD PERTURBATION EQUATIONS 

We now assume small coplanar perturbations with vi'^ = u^q' + w^^', t;^'' = v^g^ + v^l, E<'' = E[,'^ + eJ'\ = + and 
_ff'*' = J/q*^ + _ffj'' in each disc. For these perturbations in each disc, the following expressions of Fourier decomposition are used 

= (r) exp[i(a;i - raO)] , (47) 
«w = A^'' (r) exp[j(Ljt - mO)] , (48) 
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with m as an integer for azimuthal variations and i — s, g for the stellar and magnetized gaseous discs respectively; parameter uj is the 
angular frequency of coplanar perturbations and is connected to the angular wave pattern speed fip by u = mfip . The perturbation surface 
mass density is set to the following form 



(49) 



with f3i = Pi — iv 12 . The perturbation coefficient 5} is a small amplitude coefficient for each disc. The complex phase factors exp(jm^) 
and exp(ii' In r) represent azimuthal and radial variations, respectively. In general, index /3i can be different from /3 of the equilibrium disc 
configuration. The dispersion relation for coplanar perturbations in a gravitationally coupled disc configuration is then given by 

[(a(^))2 - GEi=VKn(/3?) - ojf'VVQ*''] [e(a(f')2 - eGE<«'ry„(/3f ) - aj'^VVQ^"'] = [GE<'Vy^(/3?)] [GE<^Vy™(/3?)] , (50) 



with the notation abbreviation iu'-'^ = uo ~ mQ,Q , and Q'"' and Q'®' are two abbreviations defined in Appendix|A] Derivation details of 
this dispersion relation are presented in Appendix IaI Dispersion relation dSOt is of great importance and is very powerful as it contains all 
useful information about disc dynamics in the presence of coplanar perturbations. The regime of small uj corresponds to the quasi-stationary 
situations for the QSSS hypothesis (e.g. Bertin & Lin 1996 and extensive references therein). To be specific and for simplicity, we set lj = 
in our frame of reference to study the case of non-axisymmetric stationary perturbations with m 7^ which can be either aligned or unaligned. 
Dispersion relation l |50t can be reduced to a stationary dispersion relation as 



GSq Ym {Pi 
[„{»)] 2„2^i 



D 



(s)l2 



- 2(1 - /J) 



eGSj'^Y^iPf) 
e[a(9)]V2'3 



m2 - 4(/3S*)2 + 2f3f + 213 _ 
(,).2r m^- 2(1-/3) 



D 



m? - 4(/3f)2 + 2(3f + 2/3 



r2/3 



e[a(9)]V2/3 



(51) 



where D^"'^ and are respectively the two disc rotational Mach numbers as defined in Appendix IaI 

By using equations i42t . i43t , iA37t and l lA38t , stationary dispersion relation | |5U can be cast into the explicit form of 



(l + 2/3)G(5'^«')^/<^+^^>F^(/3i^; 
(l + 4/3)fc(^)5^/(^+^« 



g(gto))l/(l+2;3)2^y^(^)(^ + go)(l + /) 



(l + 4/3)fcW5o 



1/(1+2/3) 



(2/3 + l)[m^- 2(1-/3)] 
m2 - 4(/3f)2 + 2/3f + 2/3 



(l + 2/3)6G(g^''')^/<^+"^'y„(/3f) 

(i + 4/3)efc(9) 



2/3G(5'^^')^/<^+^^)yo(/3)(/(l + 5-') + + e) 
(1 +4/3)efc(s) 



(2/3 + 1) [m^ - 2(1 - /3)] ] _ (l + 2/3)G(5'<«')i/(i+2«y„(/3f) (1 + 2/3)G(S','/>)i/(i+2«y„(/3f ) 
m2 - 4(/3f )2 + 2/3f + 2/3 



(l + 4/3)fc(=)5y(i+^« 



(1 + 4/3)efc(9) 



(52) 



Dimensionless parameters /3, m and /3i can be chosen within the allowed ranges. Dimensionless parameters e and Q depend on the isopedic 
magnetic field characterized by the dimensionless isopedic parameter A. For parameters S'^\ fc'°', fc'^' and So, typical values of late-type 
spiral galaxies are adopted and described in details in Section|6]below. 



5.1 Relation between Dark Matter Halo and Isopedic Magnetic Field 

The goal of our investigation is to explore the functional relation /(A) between the amount of dark matter and the isopedic magnetic field 
with other parameters specified. For this purpose, we further introduce several simplifying abbreviations for the derivation of /(A) relation, 
viz. 

M= " ^ ^) ^ GSj'' (53) 

m2 - 4(/3^ )2 + 2/3^ + 2/3 ' [a(^)]\^P ' 9 [ats)] \^f> ' 

which are all dimensionless. The simplified form of stationary dispersion relation J5U then appears as 

{1 - B^'^Y^ipi) - [Z)(^)]'a/} {1 - eB(«'y„(/3f) - [d'^'^]^ m] = [B^'~>Y^{l3i)\ [B'^'^Y^iPt)] , (54) 
where the two disc rotational Mach numbers squared are defined by 

[d<='] ' = 2/3B(''' Yo(/3)(l + 5o)(l + /) - 2/3 - 1 , (55) 

= 2/3B('''yo(/3) [/(I + 5,7^) + 5^^ + - 2/3 - 1 . (56) 

The two parameters M and Ym(/3i ) are complex in general as parameter /3i = /3i — ii'/2 is complex in general. But for two separate special 
cases Pi = /3i (i.e. u = Q) and /3i = 1/4 with / 0, we have shown in Appendix dAU that both M and Ym {Pi) become real numbers. 
In the literature, the situation of = is called the aligned case because perturbation patterns appear aligned along the radial direction; that 
is, only azimuthal variations but no radial variations are involved, and the situation of 7^ is called unaligned case because perturbation 
patterns appear spiral-like. 

^From now on, we shall only focus on these two special real cases 1/ — with pf = Pi — P and v ^ Q with Pi — 1/4 in the remainder 
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of this manuscript. Unless otherwise stated, we have either fif — (5 or f3f — 1/ A — iu/2 in all the following equations. With our notations, 
stationary dispersion relation l l54t for non-axisymmetric stationary perturbations can be cast into the following form of 

[1 - B^'^YmiPt) + (2/3 + 1)M - 2/3B''''Fo(/3)(1 + (5o)(l + f)M] (57) 

x{ - 2/3B(«'Fo(/3)[/(1 + 5(7') + + A^I + 1 - eB^^'^Y^Wt) + (2/3 + 1)M} = b'^'B^^' [i^n (/??)]' ■ 

In the absence of the gravitational coupling on the right-hand side (RHS) of equation J57t . the first and second factors on the left-hand side 
(LHS) represent stationary dispersion relations in stellar and magnetized gaseous discs separately. That is, for a single stellar disc, we can 
determine a / value for stationary perturbations. Or, for a single magnetized gaseous disc, we can also determine a / value for stationary 
MHD perturbations. When the two discs are coupled by gravity, we need to determine the / parameter in a joint manner by dispersion 
relation l l57t for stationary perturbation patterns in both discs with different angular rotation speeds. Again, two notational abbreviations are 
introduced below for the convenience of further derivations, viz. 

Ci = l~ 2/3B'^Vo(/3)(1 + 5o)M - B^'^Y^i/St) + (2/3 + 1)M , (58) 
C2 = l- 2f3B^'^Yo{f3){e + 5o'-)M - eB^^'^YmiPi) + (2/3 + 1)M . (59) 
Then stationary dispersion relation l |57t for coplanar perturbations can be rearranged into the form of 

[Ci - 2/3B(^'yo(/3)(l + 5o)Mf] [Ca - 2/3B(«Vo(/3)/(1 + 5o^)M] - B^"^ b'^o^ [V,n(/3?)f = . (60) 
This leads to an explicit quadratic equation in terms of the gravitational potential ratio / given below 

.2 _ CiB(^>(l + 5o-^) + C2i?<-'(l + 5o) C1C2 - [Yn^iPi)]' ^ 

■^2pB(')B(9)YoW)M{l + 5o)il + 5-^) 4(3^B(')B(3)[Yo{(3)]^]VPil + 5o){l + S^^) 

By further introducing two notational abbreviations 

y(9) = 5{9)(l + 5-i) ^ V^"'! = B^'^l + So) , (62) 

the two roots /i and /2 for the gravitational potential ratio / are simply 

Ci\/'«) + C2V-<^' ^ {[Cil/(«)-C2V^'^'f + 4[y„(/3f)]VWT/(«)i?(-)i?(^)}'/' 
■'^'^ ~ ■'^ ^ 'if3VMV(3)Yo{(3)M 4/3y(=)V(f)yo(/3)M 

where the two roots /i = /+ and /2 = /_ correspond to the plus and minus signs, respectively. All variables are real numbers in our model 
study here for either /3i = /3 or /3i = 1/4 — iiy/2, and there always exist two roots for / which can be either both real or a pair of complex 
conjugates. By our definition of / = <!>/$, we should identify the real positive solution(s) for / ratio. This leads to the requirement that the 
determinant under the square root in expression l l63t must be positive which turns out to be automatically satisfied. Another requirement is 
that in equations JA3U and (lA32} of Appendix A, both v^q and t;|o should be non-negative which leads to the following conditions, 

_ (l + 4/3)fe(-)K°T^^^'+^^^ ^ _ e(l+4/3)fc(^)[g^^)]^^/^'+^^> , (l-.)g(^) ^ 

2/3GFo(/3)5u '''''^ 2f5GYo{f3)So ^ So ^ ' 

Together with inequality J46t . there are thus total four necessary conditions on gravitational potential ratio /, viz. 

• a) / e K+ , 

. b) S^«)(/)>0, 

•c) [v%^]\f)>0, 

.d) [vitA\f)>0, 

in our model analysis for global stationary perturbation patterns. 



6 PARAMETERS FOR MAGNETIZED SPIRAL GALAXIES 

In dispersion relation l l52t . dimensionless parameters m, /3, jSf, e and B should be specified to characterize stationary perturbations in 
composite disc system. Physically, integer m indicates the number of spiral arms for a perturbation pattern, /3 is the power-law exponent 
to characterize radial variations of the different unperturbed background variables, jHf is the exponent to characterize radial variations of 
perturbations. Parameters e and B are the reduction and enhancement factors associated with an isopedic magnetic field, respectively. For 
relevant coefficient Sq'\ we adopt parameter ranges listed in Table[T]in reference to estimates from observations (e.g. Roberts 1962). Two 
parameters fc'"' and fc(9) are functions of /3 parameter and closely relate to the sound speeds (as given in Table [T) according to equations 
( l42t and ([43j with cgs unit of cm^'V(s^ g""')- Typical values for k^"^ and fc'"' are of the order ~ lO"' cm^'V(s^ g"~')- The range for 
6 is a function of (3 and So and cannot be given here simply but can be computed in a straightforward manner. For the Gamma functions 
r(z) involved in our computations, we use the code developed by Zhang & Jin (1996)0 Our study explores behaviours of potential ratio / 
dependence on isopedic parameter A for globally aligned and unaligned logarithmic stationary perturbation patterns. Physically, this study 



^ For Gamma functions r(z) in C/C++ languages of real and complex arguments, the reader is referred to website [http://www.crbond.com/math.html . 
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Table 1. Parameters for the study of late-type disc galaxies with isopedic magnetic field. Numerical ranges for 5q and Sq are calculated for /3 = [0, 0.5]; 
the ranges for A and e are determined with the given ranges of and tabulated below. 



Variables 


Value 


^0 


2 X 10""' - 2 X 10-3 g cm-2 




lO"'^ g cm--^ 




(1019, 10«)gcm2'5-i 




0.05 


a(9)(5kpc) 


1 X 10^ cm 


a(''){5kpc)^ 


3 X 10"^ cm 




(1, 10) AtG 


Bz in central regions 


(20, 40) fiG 


A 


(0.03, 3.25) 


e 


(-1110, 1) 



50 



40 



30 



20 



10 



0.5 



m=1 

m=2 



1.5 
X 



2.5 



50 
40 
30 
20 
10 



p=0-1.5o=005 



0.5 



m=1 
m=2 



1.5 

X 



2.5 



Figure 1. The dependence of the gravitational potential ratio / on the dimensionless ratio A = 2ttG ' /Bz and m in the stationary aligned case for 

late-type disc galaxies with So = 0.05. Different /3 values are studied for the aligned case of = and /3i = /3. For this and all the following figures, values 
of all other parameters, which are not explicitly mentioned, are listed in Table[T]of Section|6] For f} = 0.1 and m = 1 (right panel), / — »~ 0.2, while for 
P = 0.3 and m = 1 (left panel), / 0.5. 



reveals possible relationships of an isopedic magnetic field and the gravitational potential of an axisymmetric dark matter halo for stationary 
global perturbations. Since Vm(/?i) ~ Y-rniPi) (see footnot^ and M{vn?) — Af ((— m)^), we focus on the case of m > without any 
loss of generality. We present only the cases for m — 1 and m. = 2 since a previous analysis of Lou & Wu (2005) has shown that all curves 
for m ^ 2 bear the similar form, whereas the case of m = 1 carries a unique form. 



7 GLOBAL STATIONARY ALIGNED PERTURBATION PATTERNS WITH u = AND Pi = (3 

For global stationary aligned cases with i/ — 0, v/e have further taken a special case of /3i = /3i = /3 for perturbations carrying the same 
scale-free index as that of the background equilibrium disc configuration. Relations /(A) for different combinations of (3 and 5o are shown 
in Figures [T] and |2] as examples. We note first that no sensible /2 root can be found in Figure[T]satisfying requirements c) and d) listed after 
inequalities l l64t . In the following, we only show stationary perturbation solutions with / satisfying all four requirements a)— d) above. 

We have explored P values within the range /3 = [0.01, 0.49]. Figure[T]shows examples for two different /3 values. The numerical results 
for P > values that we have studied do not show remarkable differences among all m values. While m 2 give monotonically decreasing 
functions /(A) for all P > 0, m — 1 case leads to an almost constant /(A) whose value decreases with decreasing p. This latter trend of 
decrease cannot be easily discerned in the figures shown, but can be readily found by checking the specific numbers. For m ^ 2, stationary 
solutions for aligned perturbations can be found in all figures for small A corresponding to stronger magnetic fields. For m — 1, stationary 
solutions for aligned perturbations exist in a wide range of A values. In general, / = 10 is a typical value for spiral galaxies containing dark 



* A mathematical proof of this relation can be found in Appendix A of Wu & Lou (2006). 
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Figure 2. The aligned case with i/ = and /3i = /3 is studied. The dependence of the dark matter amount / on the dimensionless ratio A and m in the 
stationary case for a spiral galaxy with /3 = 0.3 and different So con'esponding to different epochs of a spiral galaxy evolution. 



matter halos. In Figure[T]of a late-type galaxy, / = 10 leads to A = 0.1. For a gaseous surface mass density of E*,''' = 2 x 10"" R cm , 
one gets = 4 fiG which is a realistic number and for eJ,®' — 2x 10~'^ g cm~^, the magnetic field strength is Bz = 40 fiG which may be 
reached in the central regions of spiral galaxies. It is then realistic to have stationary aligned m = 1 perturbations corresponding to lopsided 
patterns. 

An increase of So in Figure [2] should be understood as going back to earlier evolution phases of a spiral galaxy. This is because at early 
times, the gaseous disc contains more mass and during the course of evolution, the stellar disc becomes more and more massive as a result of 
star formation. When So is increased, the stationary range for m ^ 2 is also increased to larger A values whereas that of m = 1 is decreased. 
For larger values of So, m ^ 2 gives stationary solutions in the entire range of A. In addition, a second solution of m = 2 can be also found in 
the range around A G (0.7 ~ 1) for the two larger So values. In the limit of So — > oo, the stellar disc contains a negligible amount of mass as 
compared to that of the magnetized gaseous disc. Effectively, this situation can be understood as the presence of only one single isopedically 
magnetized gaseous disc which was investigated earlier by Wu & Lou (2006). For So > 0.05, we study again the case of f — 10. Here, A 
can be much larger than only 0.1 for late-type spiral galaxies. For So = 1000, one gets A — 0.4. In this case, magnetic fields are B^ = l/iG 
for eJ,"' = 2x 10"*g cm"^ which is also reahstic and B^ = lO^G for eJ,"' = 2x lO^^g cm"^. 

The case of m = 1 represents an exceptional case in our exploration so fan Observational evidence for the existence of such so-called 
lopsided galaxies can be found in Baldwin, Lynden-Bell & Sancisi (1980). For a late-type spiral galaxy with So ~ 0.05 and a positive (3 
value, the stationary solutions in Figure 1 can be roughly divided into two ranges of A. On the left side with small A or strong magnetic field 
strengths, stationary solutions for m ^ 2 are found. In almost the entire range of 0.2 ^ A ^ 3.25, the m — 1 case has stationary solutions. 
In general, in the range of weak magnetic fields only m — 1 has stationary solutions. For strong magnetic fields, only m ^ 2 can possess 
stationary solutions. The limit of A ^ oo can be regarded as the absence of magnetic fields which was analyzed by Shen & Lou (2004). 

For larger So corresponding to an earlier phase of an evolving spiral galaxy as shown in Figure 2, the stationary range of A for m = 1 
shrinks whereas that of m ^ 2 increases to all possible physical values of A; within a certain range of A, a second solution for m ^ 2 also 
comes into existence. 

Independent of parameter variations, /(A) remains always almost constant for large A or weak magnetic field strengths. This is sensible 
as weak magnetic fields are not expected to affect the entire perturbation configuration significantly (e.g. Wentzel 1963). For strong magnetic 
fields, we need to have a larger amount of dark matter to maintain stationary perturbation configurations. In the aligned case, this increase 
of / can reach values of several tens, e.g. in the limit of A = 0.03, we have / ~ 45, 55, 75 for /3 = 0.49, 0.3, 0.1, respectively. The 
exploration of different So values has revealed additional information: in disc evolution in the early Universe, magnetic field should have 
played a more important role in a spiral galaxy. This is because the range of A with a significant change of / appears larger; thus weaker 
magnetic fields should have also affected / values. Note that the magnetic field cannot cause a growth of dark matter halo in a spiral galaxy. 
Nevertheless, this analysis reveals the relation between the dark matter (/ ratio) and the isopedic magnetic field strength (A parameter). 
As expected for weak magnetic fields, one may ignore B^ in general. While for strong magnetic fields, there is only one possible / for 
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Figure 3. The function /(A, m) for (3 = 0.3, <5o = 0.05 and /3i = 1/4 in the unaligned logarithmic case with /3i = 1/4. In the lower two panels, the two 
branches of /i and /2 with /i > /2 are shown for each m ^ 2 (see Appendix B). Note that for small A, /i and /2 branches do not cross each other, although 
they appear to be extremely close. 



each A in order to maintain globally stationary perturbation configurations. Physically, a larger / ratio also leads to a larger disc rotation 
speed according to equations JAB It and iA32t . For example, one estimates u^q' = 270 km s^^ for / = 20, r = 10 kpc, /3 — 0.1, 
a'"' — 30 km s^^, Eq*' — 10^^ g cm^^ and So = 0.05 for a late-type spiral galaxy; while for / — 30 with other parameters the same, one 
obtains Wgp' = 330 km s^^. Although very large disc rotational velocities have not been observed, the part with large / can be regarded as 
theoretically plausible solutions for spiral galaxies. 



8 GLOBAL STATIONARY UNALIGNED LOGARITHMIC SPIRAL CONFIGURATIONS WITH v^Q AND /3i = 1/4 

When MHD density wave perturbations propagate in both radial and azimuthal directions, we have global unaligned logarithmic spiral 
patterns with v ^ Q and m 7^ 0. In our consideration here, we take Px — I3i — iv/2 7^ /3 for 7^ 0. In this case, M and Y'o(/3i) as 
defined by equations l l53t and l l34t are complex in general, while dimensionless parameters B'^"^ and S'*-* as defined by equations l l53t are 
always real. In the special case of /3i = 1/4, we have shown that both M and Yo{f3f) are real (see Appendix A). Physically, this case of 
/3i = 1/4 corresponds to a constant radial flux of angular momentum (e.g. Goldreich & Tremaine 1979). Again, the stationary condition 
for logarithmic spiral patterns implies a certain relation between / ratio (related to the dark matter halo) and A parameter (related to the 
isopedic magnetic field). For i/ 7^ 0, a special attention is further paid to i/— dependence of /(A) (see footnote 3). For different 1/ values 
(related to the radial wave number) in a warm disc, the curves /(A) are shown in Figure |3] For small v, at least one stationary solution 
exists for each m > 0. However for increasing v towards the tight-winding regime, stationary perturbation solutions of m = 1 are shifted 
towards weaker magnetic fields. For v — 50, solutions for m = 1 are found within the range A ~ (1.1, 3.25) and for u = 100, they are 
within A « (1.2, 3.25). Meanwhile, stationary perturbation solutions for m ^ 2 grow into a larger range of A. If v is increased further 
(i.e. towards the extremely tight-winding regime), two stationary perturbation solutions exist corresponding to each m. Thus according to 
expression l |63t for the two roots of /, /i and /2 are the upper and lower curves in Figure |3] respectively. In the range of A for stationary 
perturbation solutions, /i remains more or less constant so that a change of magnetic field strength appears independent of the amount of 
dark matter. Meanwhile, /2 changes over a larger range of A, indicating that the magnetic field and /2 relate each other more closely. The 
magnetic field is directly connected to A. If /i remains constant for all A, then the reverse direction of reasoning cannot be used, i.e., for a 
given /i — 19 and u — 50, the exact value of A and so Bz cannot be determined with certainty. In contrast to /i, /2 varies for = 50 over 
a larger range of A = (0.1, 0.5). In this range of A, one can determine for each /2 the associated A and thus B^ through B^ = 2tv\^T,q/ X. 
The tighter the spiral arms are wound, the more dark matter is needed in order to maintain stationary logarithmic spiral patterns. This is a 
general trend of variation and is valid for both roots of /. Nevertheless, there exists a certain i/q value such that the second solution /2 appears 
with a much less amount of dark matter for u > vq. In other words, with increasing f the dark matter amount also increases, except for the 
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second solution with v > vq. Physically, the two roots /i and /2 are equally valid. But compared with the observational value of / ~ 10 
for typical spiral galaxies, /i appears too large, e.g. for = 50 and /i ~ 19 within the range of A = (0.2, 3.25). In comparison, the lower 
f-z appears more plausible because of its smaller values instead of the very large values of /i. For v = 100, one can again find a relation 
between the fz curve and the magnetic field (represented by A) and it is thus also possible to determine Bz through the / ratio. If the curve 
/2 (A) is sufficiently steep, so that a one-to-one relation is present between fz and A, then the magnetic field strength can be determined 
by integral l |18l l into 

B. ^ . (65) 

The only additional variable that one needs to know is the surface mass density of the magnetized gas disc. In this paper, we have not explored 
situations of v and /3i 7^ 1/4 for which both parameters M and Fm(/3i — if/ 2) are complex. In this case, an additional study is needed 
to explore the variation of /3i and its influence on the gravitational potential ratio /. Physically, /3i is the constant scale-free index of surface 
mass density perturbations. A larger Pi corresponds to a steeper radial profile of stationary coplanar perturbations in magnitudes. 



9 MODEL RESULTS FOR TWO LIMITING CASES 

In the first limit of large A regime (i.e. Bz —> 0), we readily reproduce the model results of Shen & Lou (2004) where the gravitational 
coupling of two scale-free stellar and gaseous discs in the absence of magnetic field is analyzed. For the absence of magnetic field, we choose 
a sufficiently large isopedic parameter A = 24. The dependence of ratio / on various parameters is shown in Figure |4] For the choice of 
So = 0.05 and v — 0, only the case of m = 1 has stationary perturbation solutions. The variation of So enables m ^ 2 to have stationary 
perturbation configurations as already mentioned in Section|7] For 1^ = 0, each m corresponds to only one stationary perturbation solution. 
When V is sufficiently small, only m — 1 perturbation configurations have stationary solutions. As 1/ is gradually increased, one can see 
clearly in Fig.|4]panel c) that stationary perturbation solutions for m ^ 2 also appear. 

In the second limit of large So (i.e. Sq*' 0) for young galaxies in the early Universe, we recover the results of Lou & Wu (2005) where 
global MHD perturbation configurations in a single isopedically magnetized gaseous disc are studied. Figure|5]shows different variations of 
/ ratio for the aligned and unaligned logarithmic cases in that regime. In principle, all m values can have corresponding stationary solutions 
for global perturbation configurations. This depends on the choice of A and the combination of different parameters such as k^"\ fc'"' or 5'^' 
and so forth. For m = 1, the starting point of stationary solutions for the aligned case is A ~ 1.1 and for the unaligned case it is slightly 
shifted to A « 1.2. For m — 1 the range of ratio / = (0.01, 0.5) is very small for v — sls compared to range / = (1, 60) for f = 100. 
For m ^ 2, a variation of /3 parameter in the range of f3 — (0.1, 0.49) for the background composite disc configuration does not affect 
features of /(A) curves very much. In the model of Wu & Lou (2006), potential ratio / is studied for different properties of the magnetized 
gaseous disc alone (such as the sound speed) at a certain specified value of A parameter. The special case of one single magnetized gaseous 
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Figure 5. The limiting case of only one single isopedically magnetized gaseous disc is approximately achieved by setting Sq = 1000. Two branches f\ and 
/2 roots are shown with /i > /2 for each m ^ 2 (see Appendix B). While they appear very close, /i curve does not intersect /2 curve at small A. 



disc in our model is complementary to their analysis. Again, the basic fact that typically / < 10 must be kept in mind for applications to 
spiral galaxies. Due to their smaller values for v = 100, the minus-sign solution /2 is regarded as more plausible. 



10 CONCLUSIONS AND DISCUSSION 

We have investigated global stationary solutions for aligned and unaligned logarithmic perturbation configurations with a constant radial 
flux of angular momentum and paid special attention to the roles of dark matter halo represented by the / parameter for the gravitational 
potential ratio and the isopedic parameter A. For this purpose, the stationary dispersion relation for global perturbations is derived and a 
quadratic equation of / is obtained. In our model formulation, the stationary assumption is applied as a very special limiting case of the 
QSSS hypothesis (e.g. Bertin & Lin 1996 and references therein) which itself is already a very strong requirement to the disc dynamics. For 
a non-vanishing small pattern speed of a few km s^^kpc^^ (adopted in the QSSS theory), a similar relation between / and may be also 
derivedlfl Our sample calculations show that there are three possibilities for potential ratio /, corresponding to no solution, one solution and 
two solutions, respectively. In general, potential ratio / is a fairly complicated function of several independent model parameters involved. In 
order to explore the dependence of potential ratio / upon /3, Sq, \,m and u, we have chosen typical values for these relevant parameters to 
characterize late-type disc galaxies. Therefore, these results should be applicable to typical spiral galaxies for comparison. The main results 
are now summarized below. 

10.1 Global Aligned Perturbation Configurations with v — Q and f3i = f3 

Global aligned perturbation configurations with 1^ = are studied for typical late-type galaxies with different /3 values, corresponding 
to various power-law radial fall-offs of the disc rotation curve. Examples of our numerical exploration show that in the estimated range 
of A for disc galaxies, a variation of /? parameter has no significant impact on global stationary perturbation configurations with m ^ 1. 
As a general trend, a stronger magnetic field allows existence of stationary perturbation configurations for rn ^ 2. For weaker magnetic 
fields, only stationary perturbation configurations with m — 1 (i.e. lopsided cases) may exist. We have also explored stationary perturbation 
configurations for a disc galaxy at different epochs of evolution by varying the disc surface mass density ratio Sq = ®' /Eq*' . By our sample 
calculations, the variation of So strongly affects the behaviour of /(A). This is a consequence of the fact that the magnetic field directly affects 
the gaseous disc and in an earlier stage of galaxy evolution, the gas fraction is higher in a composite disc system. With increasing So ratio 
towards earlier epochs ( i.e. at higher cosmological redshifts in the early Universe), the range of allowed stationary perturbation configurations 
for m ^ 2 also increases while that for m — 1 shrinks. In our simple scenario. So parameter marks the evolution of a disc galaxy in the 
expanding Universe. At the beginning, a proto-spiral galaxy is presumably composed of a nebulous gas disc and a dark matter halo (e.g. 
Lou & Wu 2005; Wu & Lou 2006). During the evolution of a spiral galaxy, more and more stars are born and die and thus the disc mass 
ratio So decreases slowly with time. With this scenario in mind, we interpret Figure|2]as follows. At an earlier stage of disc galaxy evolution, 
stationary solutions of global perturbation configurations with m ^ 2 exist over a wider range of A. The influence of magnetic fields there is 
also stronger as the curve /(A) there varies significantly over a larger range of A, indicating a wider range for magnetic field strengths. 

This brings out an extremely interesting evolutionary perspective for speculations. Observationally, the global star formation rate (SFR) 
in spiral galaxies is an important indicative parameter to characterize the galactic evolution. Conceptually, if one attempts to relate this 
SFR with disc instabilities, the application of Toomre's criterion (Safonov 1960; Toomre 1964) to a stellar disc alone would be insufficient 
because stars are directly born in the gaseous disc. Following this line of reasoning, one needs at least to explore instabilities in a composite 

^ Wu & Lou (2009 in preparation) obtained a dispersion relation for non-zero but small pattern speed and in the absence of the stellar disc component (see 
also our expression 150) and Wu & Lou 2006 for very small angular perturbation frequency uj). The range of potential ratio / should be determined in order to 
evaluate the importance of the isopedic magnetic field. Intuitively, the results for quasi-stationary and stationary configurations are expected to be qualitatively 
similar. Quantitative deviations are expected to be proportional to certain powers of pattern speed (e.g. Wu & Lou 2006). 
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system of stellar and gas discs (e.g. Lou & Fan 1998, 2000 and references therein) in the presence of a massive dark matter halo. The 
importance of magnetic fields is generally recognized in the dynamics of star-forming cloud cores. Nevertheless, how to specifically relate 
physical processes of star formation in clouds on much smaller scales and of large-scale disc instabilities remains a challenging problem 
due to the tremendous differences in scales. In spite of this challenge, we have developed intuitive feelings that regions of high gas density 
and strong magnetic fields on large scales are expected to be vulnerable or favorable to active star formation processes. Now in a highly 
simplified dynamic manner, our model analysis brings together several important aspects for this physical consideration, viz. dark matter 
halo, magnetic field, and higher gas fraction (i.e. larger 5q) in a composite disc system. The initial conditions for forming a proto-galactic 
disc nebula such as the dark matter halo, gas disc and magnetic field are expected to be statistical with fluctuations. The instability properties 
of such a magnetized disc system in the presence of a dark matter halo may then grossly determine the global SFR and thus the galactic 
evolution. In other words, different initial conditions can lead to different kinds of evolutionary tracks. In particular, it is conceivable that 
the SFR may have highs and lows along galactic evolutionary paths in the expanding Universe. This information would be valuable for 
understanding the overall cosmological evolution in terms of global SFRs of spiral galaxies. 

To compare a late-type spiral galaxy with an early-type spiral galaxy, we emphasize that stationary perturbation solutions for m ^ 2 
are very different in these galaxies. In a late-type spiral galaxy of So ~ 0.05, stationary perturbation configurations of m ^ 2 can only exist 
for strong magnetic fields while in an earlier stage of evolution, such solutions can exist for almost all A values. These results suggest that 
multiple-armed disc galaxies may be more numerous in the early Universe. It may be possible for disc galaxies to change patterns during 
the course of their evolution. Conceivably evolving on cosmological timescales, a disc galaxy pattern may alternately become stationary or 
quasi-stationary as evidenced by the fact of changing So and A, e.g. / and m ^ 2 are shifted towards smaller A for stationary configurations. 
This offers a novel perspective to relate global SFRs of disc galaxies and the galactic pattern speed evolution. For example in numerical 
simulations, one may start from initially stationary perturbation configurations and then explore time-dependent quasi-stationary behaviours 
as well as nonlinear effects of a composite disc system by adjusting one or several relevant parameters systematically. 

10.2 Global Stationary Unaligned Logarithmic Spiral Conflgurations with /3i = 1/4 and v ^ Q 

Parameter characterizing radial variations of coplanar perturbations affects / ratio in many ways. For example, in spiral galaxies with 
more tightly wound logarithmic spiral arms (i.e. larger v), more dark matter is needed in the massive halo to sustain stationary perturbation 
configurations. Furthermore, with increasing u towards the WKBJ regime, each / may correspond to multiple stationary configurations. 
For very small v in the opposite limit, only perturbation configurations with m = 1 have stationary solutions (i.e. lopsided configurations). 
Stationary perturbations with m = 1 also have solutions for larger v, but only in a certain range of A ~ (1-2, 3.25). Larger values have 
one stationary solution for all m ^ 2 (see Appendix B), namely the /i root of expression l l63t . If v is increased further, then the second root 
/2 of expression l |63t also fulfills all necessary physical requirements on / parameter. Here, /2 root is always smaller than /i root and we 
expect that whenever there are two theoretical solutions possible, then /2 tends to be the more realistic root because of its relatively lower 
values. The limiting case of — > cx) represents tightly wound spiral arms and the well-known WKBJ approximation (Lin & Shu 1964, 1966) 
should become applicable. 

The curve /2 (A) shows that there is a one-to-one correspondence between /2 root and strong magnetic fields. Therefore, by determining 
/2 ratio through observations (e.g. rotation curves or gravitational lensing effects) and equation l l40t . the magnetic field strength can also be 
determined. For weak magnetic fields, root /2 does not vary significantly enough and thus, isopedic magnetic field strength cannot be 
sensibly determined as a result of uncertainties in galactic observations. 

10.3 Two Limiting Parameter Regimes 

The two limiting cases of A ^ cxa and (5o — > oo have been explored in the previous section. In the limit of A ^ oo, magnetic field almost 
vanishes and one comes back to two gravitationally coupled hydrodynamic disc system with an axisymmetric dark matter halo. For late- 
type spiral galaxies, only perturbation configurations of m = 1 lead to stationary aligned solutions with a positive / ratio. For unaligned 
logarithmic perturbation configurations with m ^ 2, the existence of stationary solutions also depends upon the choice of dimensionless 
'radial wavenumber' v. For early-type spiral galaxies of higher 5o, global stationary solutions of all m for perturbation configurations can be 
found. As already noted, this implies the possibility of more numerous multiple-armed spiral galaxies in the early Universe. 

For So ^ oo (see Fig. 5 for large So), stationary perturbation solutions for m ^ 2 can be found for every chosen parameter set. For 
m = 1, solutions exist in a certain range of A. For example for v = 100 and So ~ 1000, the case of ?n = 1 has / roots for A — (1.2, 3.25) 
(see the right panel of Fig. 5). 

To conclude for So oo, for small v or v = Q, all m values have one / root. But for larger v, the case of m ^ 2 can have two / roots. 
This is the same as that of So = 0.05. The difference lies in the values of /. The larger the 'radial wavenumber' v value is, the larger are the 
potential ratios /i and /2 . 

For late-type spiral galaxies and unaligned logarithmic cases, a weak magnetic field does not play a significant role as expected. The 
younger the spiral galaxy is, the more important the role of a magnetic field especially in the regime of stronger magnetic fields. For very 
strong magnetic fields, it appears that / ratio is more closely connected to the magnetic field and the effect of magnetic fields cannot be 
ignored. But for weak magnetic fields, / ratio remains more or less constant for all types of spiral galaxies. 

Our results can be summarized as follows. First of all, we emphasize that the lopsided global configuration m = 1 is an exceptional 
case for which the curve /(A) has a distinctly different shape as compared to those for the cases m ^ 2 (see Appendix B). For both aligned 
and unaligned logarithmic spiral cases with m ^ 2, strong magnetic fields (still realistic in spiral galaxies) can bear a significant relation to 
the dark matter halo in maintaining globally stationary perturbation configurations. 
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Physically, we conclude that in spiral galaxies with or without radial variations in perturbations, weak magnetic fields do not influence 
stationary perturbation configurations. But when the magnetic field strength is increased, a globally stationary perturbation configuration 
would then become non-stationary if the / ratio is not high enough for the increased magnetic field strength. A change of perturbation patterns 
might be possible during the course of the galaxy evolution. For example, flocculent galaxies might represent the transitional phase for global 
pattern changes of galaxies. This prediction may be tested by numerical simulations and by deep survey of morphological observations for 
galaxies in the early Universe. 

Our model analysis has shown that the magnetic field needs to be sensibly chosen with other given disc parameters in order to maintain 
a global perturbation configuration with a stationary pattern in our frame of reference. Therefore, a variation or an adjustment of Bz may 
put a stationary perturbation configuration into a non-stationary one or a non-stationary perturbation configuration into a stationary one 
because there exists a one-to-one correspondence between the potential ratio / and a sufficiently strong magnetic field Bz for globally 
stationary perturbation configurations. For such a stationary disc balance in general, the potential ratio / can vary in a considerably large 
range depending on the choice of other relevant parameters in sensible regimes. As expected on intuitive ground, sufficiently weak magnetic 
fields exert fairly small influence on stationary perturbation configurations in our composite model. The main reason is that the magnetic 
pressure and tension together are not strong enough as compared with other forces which play important dynamic roles in our composite 
system. 

In the regime of weak magnetic fields, there also exist instabilities widely explored in the literature. The magneto-Jeans instability 
(MJI) is an instability which is based on background in-plane magnetic fields. In recent years, several authors (e.g. Kim & Ostriker 2001; 
Kim, Ostriker & Stone 2002; Shetty & Ostriker 2006) have performed numerical MHD simulations for galaxies. In particular, they studied 
azimuthal magnetic fields that lie in the disc plane (see also Lou & Zou 2004, 2006; Lou & Bai 2006). For our study of isopedic magnetic 
fields, as already mentioned, Lou & Wu (2005) have proved that a constant A is a consequence of the frozen-in condition on magnetic flux (see 
also Wu & Lou 2006). Our composite MHD model offers a two-dimensional description of an isopedic magnetic field. Due to the idealization 
of razor-thin discs, perturbations lie in the galactic plane and have no 2— components. Hence, a treatment of the magnetorotational instability 
(MRI) is impossible due to its necessary requirement of a perturbation e''^^^ along 2— direction (Chandrasekhar 1960; Balbus & Hawley 
1991, 1998; Balbus 2003). When switching over to discs of finite thicknesses, the MRI must be included and a larger / is probably needed in 
order to counteract MRIs. With our theory, the range of strong vertical magnetic fields is well studied, whereas the effect of weak magnetic 
fields is underestimated since we do not account for MRIs which play an important role for weak magnetic fields. The results that we obtain 
here cannot be directly applied to coplanar magnetic fields, thus we can make no predictions for the MJIs. 

The galactic application of our disc model results leads to two methods which may be utilised to determine either the isopedic magnetic 
field strength B^ or the mass of an axisymmetric dark matter halo. First, by using equation ( 140b through observations of a disc galaxy, one 
can estimate the gravitational potential ratio /. In the case of a stationary pattern of a perturbation configuration, only one distribution of 
magnetic field strengths is possible for this / value. This is a new approach of determining the distribution of isopedic magnetic field strengths 
in disc galaxies. Proceeding in the opposite direction, these model results may also be applied to determine / ratio by using observationally 
inferred distribution of magnetic field strengths. This appears to be an alternative method to estimate the halo mass of dark matter. The 
real challenge of these procedures is to determine whether a perturbation pattern is stationary or quasi-stationary through independent 
observational diagnostics. In practice, the relevant parameters that need to be determined by galactic observations for our proposed method to 
work are m, /3 (the scaling index for radial variations of unperturbed disc variables), /3f (the scaling index for radial variations of perturbation 
disc variables), unperturbed disc surface mass densities Sq*' (r), effective sound speeds a'*', either A (for the determination of gravitational 
potential ratio /) or the other way around / (for the determination of distribution of magnetic field strengths Bz)- It is indeed a challenging 
task for observations to specifically identify some of these parameters and as far as we know, galactic observations so far still have not 
determined completely all these parameters with error bars in one galaxy. 

In addition, this work can also be understood as a preparation for the following study on MHD density waves in such a composite disc 
system since all relevant quantities and their relations among each other are presented in this work. In reference to singular isothermal disc 
(SID) models of Shu & Li (1997), Shu et al. (2000), Lou & Shen (2003), Lou & Zou (2004, 2006), Shen et al. (2005), there exist two classes 
of solutions for stationary magnetohydrodynamic (MHD) perturbation configurations with in-phase and out-of-phase density perturbations 
in the two discs. We expect for the case of a scale-free stellar disc and an isopedically magnetized scale-free gaseous disc embedded in an 
axisymmetric dark matter halo, that there are also two classes of in-phase and out-of-phase density perturbations. The formulae and results 
obtained here can be used later to analyze such MHD density waves. 
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APPENDIX A: DERIVATIONS OF STATIONARY DISPERSION RELATION AND RELEVANT EQUATIONS 



For coplanar MHD perturbations in a background rotating axisymmetric disc system with — v).q + v).''( 



„{') 



Va' = V, 



(0 



E^*' = S[,'' + E^*' , = + and H'-^^ = H,^'' + //f ' in equations - (l9j, where superscript i within parentheses can be 
set as s and g for stellar and magnetized gaseous discs respectively, it is straightforward to obtain linearized perturbation equations below 
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where u'o' = for the background disc configuration. Equation dAll l represents the two mass conservation equations, equations l lA2b and 
iA3l are the two radial momentum equations, and equations ( IA4b and l lASb are the two azimuthal momentum equations; effects of the 
isopedic magnetic field are subsumed in the two parameters O and e (see equationsl22land[23t for the magnetized gaseous disc in momentum 
equations iA3\ and iA5i ; and the gravitational coupling between the stellar and gaseous discs is represented by the sum of gravitational 
potential perturbations or ^j"'' +e"l?j®' in two-dimensional perturbation equations ( IA2b to JA5t . As the axisymmetric gravitational 

potential of a massive dark matter halo is presumed to be unperturbed for simplicity, its perturbation does not appear on the right-hand sides 
of equations iA2l to dASb FI For perturbations cast in the form of Fourier decomposition, we write for example 



„W = A'^'> (r) exp[i(Lot - m9)] , (A6) 

= B*'' (r) exp[i{ujt - m6)] , (A7) 

where m is an integral number to characterize azimuthal variations, lu is the angular perturbation frequency, A'*-' (r) and B'*-' (r) are small- 
magnitude functions of radius r, and superscripts i = s, g stand for associations with stellar and gaseous discs, respectively. Two new 
frequency parameters are now introduced below for simplicity. The first one is the abbreviation Zj''' = 
so-called epicyclic frequency Kq*' of a rotating disc defined by 
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Substituting all expressions of Fourier harmonics into equations JAlb — l lASl l. we readily obtain the following perturbation equations 
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for both stellar and isopedically magnetized gaseous discs coupled by gravity. Furthermore, the perturbation of surface mass density is set to 
the following form of Fourier harmonics 



eJ'> = V"'"!-' exp[i(c^i - m9 + u\nr)] = sfV"^"?"' exp[i(a;i - mO)] 



with S^^ being small-magnitude coefficients to justify the perturbation approach and 
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° The assumption of an unperturbed axisymmetric dark matter halo corresponds to a limiting case. It would be interesting to explore the case where the dark 
matter halo is perturbed and coupled to disc perturbations gravitationally. 
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being complex in general for 7^ 0. The unit complex factor exp(— im^) represents periodic azimuthal variations and the unit complex 
factor exp(ii^ In r) represents radial variations with v closely related to the radial wavenumber. In general, /3i parameter for perturbations 
can be different from 13 parameter used for characterizing the background axisymmetric equilibrium disc in rotation. According to several 
formulae derived by Qian (1992), the corresponding perturbation of gravitational potential associated with a surface mass density perturbation 
is simply given by 

= -GF„(/3?)5j'V-'''i exp[j(c^t - me)] = ^GrYm{pt)T.t^ (A16) 
analytically and the corresponding perturbation of enthalpy is given by 
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where 5\ = eJ^*' /E j"' is a constant ratio for surface mass density perturbations by expression ( IA14I >. and A''-*-' and K'^^^ defined by equations 
1 IAI8I 1 and l |A19t respectively are two dimensionless constants because and rEg'' oc r . For coplanar perturbations in the 

stellar disc in rotation, we then have the following perturbation equations 

v^^^iLo^^'^ - 2n<;'i;(^i' = -GK'-'^-^ [e'^V] , (A21) 

v'g\\zJ^''^ + vif -LJ = imGE'''' A-f^' . (A22) 
Meanwhile for coplanar MHD perturbations in the isopedically magnetized gaseous disc in rotation, we obtain in parallel 
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^From equations (|A2I|I — HA22|) and HA24|) — HA25|) . we directly derive expressions for radial and azimuthal velocity perturbations and 
as given below. 



G^« ^W^) ^^^^^ 



Substituting these expressions into the mass conservation equations and using expression JA14I (. we finally arrive at 
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for both stellar and magnetized gaseous discs, respectively, and for defining the corresponding Q^*'. A combination of these two equations 
for perturbations in stellar and gaseous discs leads to the dispersion relation of the gravity coupled disc configuration, namely 

I [a^'^f ~ GE(^Vf„(/3J') -u'^VVq'^'} {e [a(«']' - eGE{,^Vy™(/3f ) ~ aj'^VVQ'"' } = GE(='ry„(/3?)GE(«Vy,„(/3f ) . 

(A29) 

Relation ( |A29t is comprehensive and it contains all useful information about the disc dynamics after a perturbation has arisen there. Depend- 
ing on the purpose of investigation, one can define different parameter regimes by using this dispersion relation. We study non-axisymmetric 
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stationary perturbations with m 7^ and uo = Q which can be either aligned configurations or unaligned logarithmic spirals. For the stationary 
case of w = 0, equation l |A29t reduces to 
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where, according to expressions ( I40t and ( 141b for the background equilibrium disc system, we have disc angular rotation speeds given by 
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As the gravitational potential ratio / should be the same in both equations ( |40| ( and ( I41l >. a relation between [w^q'] ^ and [^"90']^ '-^'^ therefore 
be established, namely 
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Using expressions = fey'V and Eq*' = S'^^r ^ for the background scale-free disc system, equation ( |A33t can be written as 
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without radial dependence. We now introduce the rotational Mach number D'"^ for the stellar disc and the rotational magnetosonic Mach 
number D'^' for the isopedically magnetized gaseous disc, respectively. 



jj(,s) _ "eo _ "0 



D 



(9) 



,,(9) 



Using background disc equilibrium equations l l40t . | |4U and l lA34t . the following relations can then be derived. 
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with ratio A = O [a'*'] ^ / [a'°^] ^ ■ With these defined quantities and relations iA31\ . iA38t and l |A39t , stationary dispersion relation l lA30t 
can be cast into 



^ G4-V^(/3f) r^(,)V 



-2(1-/3) 



^ _ eGSg^ YmjPi) _ r£)(9)l 

e[Q,(9)]V2/3 L J 



m2 - 4(/3f )2 + 2/3f + 2/3 
- 2(1 - 13) 



m2 - 4(/3f )2 + 2/3f + 2/3 



Gg^°'y„(/30 Gg^''¥^(/3f) 
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Al Two Cases of Real Stationary Dispersion Relation and Relevant Mathematical Expressions 

Our main purpose is to establish relations between dark matter halo (represented by / ratio) and an isopedic magnetic field in gaseous disc 
(represented by the isopedic ratio A). For the two cases of pf — /3i — /3 (with ly — 0) and of /3i — 1/4 (with u 7^ 0), we show in the 
following that M and Ym as defined by equations l l53t and i35l become real numbers. In the literature, /3f = /3 is called the aligned case 
since the radial wavenumber parameter v vanishes in this case and thus, only azimuthal and no radial wave variations are present. In this 
case, we have 

+ 4/3(1 - /3) ' 
Y (a) - ^r(m/2-/3+l/2)r(m/2 + /3) 

"^"^^^ - r(m/2 ~(3 + l)r(m/2 + (3+ 1/2) ' ^""^^^ 
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For 7^ and thus a complex fif, the general expressions for M and Ym{P\) are 



M = 



m^- 2(1-/3) 



m2 + 2(/3 + /3i - 2/32 + 1.2/2) + ij.(4/3i „ 1) ' 

7rr(m/2 - /3i + iu/2 + l/2)r(m/2 + /3i - i!^/2) 



Y (B 

r(m/2 - /3i+ii'/2 + l)r(m/2 + /3i - jz//2 + 1/2) ' 
By setting Pi — 1/4 in expressions l lA43t and l lA44t . we immediately have 

- 2(1 - /3) 
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m2 + 2/3 + l/4 + !/2 ' 

r(m/2 + l/4 + i!//2) 



■^"^^^^^ r(m/2 + 3/4 + i/./2) 
Expression | |A46I > is based on the relation T{z)* 



two cases (3t = Pi = P (with v 
other relevant model parameters. 



0) and of Px 



(A43) 
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(A45) 
(A46) 



r(z*) for Gamma-function where * denotes the complex conjugate operation. For these 
= 1/4 (with 7^ in general), we explore numerically the behaviour /(A) by specifying 



APPENDIX B: DEPENDENCE OF GRAVITY POTENTIAL RATIO F ON PERTURBATION ORDERS M 

Here, we briefly describe the gravitational potential ratio / as a function for different m values and provide an upper limit for the difference 
of /(m + 2) — /(m). In general, the function / depends on m through functional parameters M and Ym{P\)- In the limit of m — > 00, we 
have M 1 and therefore the difference AM = M{m + 1) — M{m) — > 0. For the property of functional parameter Ym{Pi), we first 
discuss the property of Gamma-function. For Gamma functions, the well-known recursion formula is simply 



r(2; + i) =xr{x) 

which immediately leads to the following four relations 

r((m + 2)/2 - Pt + 1/2) = (m/2 - pf + l/2)r(m/2 - pt + 1/2) , 
r((m + 2)/2 + Pt) = (m/2 + /3?)r(m/2 + pf) , 
r((m + 2)/2 - /3? + 1) = (m/2 - pf + l)r(m/2 - /3? + 1) , 
r((m + 2)/2 + Pt + 1/2) = (m/2 + pf + l/2)r(m/2 + pt + 1/2) 

as applied to our case under consideration. Therefore, we have 

(m - 2pt + 1) (m 



Y 



2(.Pt 



- 2Pt] 



2Pf + 2) (m + 2Pf + 1) 



Yn^iPt 



(Bl) 

(B2) 
(B3) 
(B4) 
(B5) 

(B6) 



For m ^ c», it foOows that Ym+2{Pt) Ym{Pt), indicating the difference AY^iPf) = y„+2(/3?) - Ym{Pt) ^ 0. Based on 
these analyses, we conclude that for very large values of m, there is no significant difference between /(m + 2) and /(m), and therefore 
/(m + 2) — > /(m) in the limit of large m values. 
For m = 2, we have from equations l l53t and I 
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1 + P 
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Y,{pi)~Y2{Pt) 
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Y2m 



■i + pt- 2{pif 



a + pf- 2{PfY 



Y2{pt) = NY2{pt) 



(B7) 
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which defines the coefficient A'^ as a function of pf that is complex for u ^ 0. The difference Af(4) — M{2) is a function of both /3 and pf. 
The largest difference is found for P = -0.25 and pf ~ 0.2 - 0.3 with AM = A/(4) - M{2) = 0.46 (see Figure|BTJ. For applications 
to spiral galaxies, only positive P values are physically relevant. But for a theoretical consideration, a negative P is allowed to estimate limit. 
To conclude, we obtain 

AM s; 0.46 . (B9) 

For a real range of pf G [—0.5,0.5] with ly = 0, functional parameter N{Pi) is plotted in Figure lB2l The smallest A'^ is —0.6 for pf = —0.5. 
From equation ( |B6b , one can also see that Ym+2{Pi) is always smaller than Ym{Pf)- This fact explains the negative values shown in Figure 
IB2I For the absolute difference, we get inequality 



Ay I _ ly 



Yra\ 0.6y„ 



(BIO) 



In contrast to AM which can be determined as a constant, AYm is a function of Ym ■ For A/, the upper limit is /(4) — /(2) which is different 
each time depending on the chosen parameters. As a result, the difference of two functions fin) and /(n + 2) with n = m, m + 1, • • ■ goes 
towards zero, namely 
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Figure Bl. For the left panel, the difference of M(4) — M(2) is shown as a function of /3 parameter with a fixed real jSf, with a maximum at /3 = —0.25 
and (if ~ 0.25; for a real value of /3j , the difference M(4) — M{2) appears to be a monotonically decreasing function with increasing /3 value. For the right 
panel, the difference of M (4) — M{2) is shown as a function of real fif values with a fixed value /3 = —0.25. A maximum value can be readily identified. 
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Figure B2. Variation of coefficient Af(/3j ) defined by equation jB8) as a function of real /JJ* values. The factor N{f3f) in relation Y4, = (N + 1)Y2 as a 
function of pf. 



/(m + 3)-/(m+l)^0 



for 



(B12) 



APPENDIX C: PERTURBATION DENSITY RATIO 5i IN THE MODEL DISC SYSTEM 

In order to compare the mass density perturbations in the two discs, the ratio 5\ — eJ^^'/eJ^''' is derived in the following by using equation 
( |A28I >. This ratio Si is contained implicitly in Jf'*' expressions. Using this relation for the stellar disc and oj^^' — 0, we derive 



-[i3W]'[A(»)]'[m2- 2(1-/3)] 



+ ■ 



E^'') GS<='y,„(/?f)[m2-4(/3^)2 + 2/3^ + 2/3] GS'^"'F„(/3f ) 
We substitute expression l |A39| l of D'*' into equation l lClb to relate the two discs. With A'^ = — (m^ — 2 -f 2/9) , we then obtain 

^ (D(»))2e(A(s')2 + (2/3 + 1) [e(A(''')2 - [A^'^f] + 2/3GYi,(/3)(l - e)^^^' (A^'^f 



Si = N- 



GYm{f3f)S^'^ [m2 - 4{(3fy + 2f3f + 2(3] 



(CI) 
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The ratio Si — Ej^'/E^^^' is negative in the case of out-of-phase density wave perturbations. For this case, the gravity effect is weaker and 
the MHD density wave speed is faster. In the case of in-phase density wave perturbations, this ratio 5i is positive, the gravity effect is stronger 
and the MHD density wave speed is slower (Lou & Fan 1998). 



